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, Abstract. When a black hole accretes slowly, the radiative cooling of 

the infalling gas is weak and the accretion disk does not form. A hot 
^p. collisionless quasi-spherical halo is formed instead. The properties of 

such halos are discussed. The rate of accretion, the radiative efficiency, 
and the temperature and density near the hole are evaluated. 



< 



> , 

! 1. Introduction 

(N 

^ i In 1970, Victor Shvartsman introduced the idea of hot halos around slowly- 

■ accreting black holes. Today, we believe that these halos do exist around the 

Q> i supermassive black holes in the centers of non-AGN galaxies. The radio source 

On ■ Sgr A* is the nearest example. 

The radiative efficiency of collisionless halos is determined by the electron 

Qh' heating mechanisms (§2). The rate of accretion is not necessarily given by the 

Q . Bondi formula (§3). The plasma density and temperature profiles near the black 

H ' hole are estimated in S4. We conclude in S5. 
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2. The Radiative Efficiency 



At low densities, the radiation is weak and the accreting gas is unable to cool. 
a I More precisely, the protons are unable to cool. The proton temperatures are 

close to O.lGeV near the last stable orbit. For realistic densities, the plasma 
is collisionless. Electrons synchrotron cool efficiently and are much colder than 
ions. Electrons radiate away nearly all the heat they receive. The radiative 
efficiency of accretion {rj = L/ Mc^ , L is the luminosity, M is the accretion rate) 
is determined by the electron heating processes. 

The electron heating processes depend on the magnetic field strength. The 
magnetic field should be close to equipartition with the gravitational energy. 
Shvartsman (1971) explains this as follows. For a nearly spherical accretion, 
assuming perfect flux- freezing, the radial magnetic field is oc r~^, where r is the 
radius. The magnetic energy density is cx r~^. Assuming a free fall, the 
velocity is oc r"-^/^. From continuity, the density is oc r~^/^. The gravitational 
energy density is eg oc r~^/^ <C Cm- But the magnetic energy cannot exceed the 
gravitational energy, therefore Cfyi ~ Cq. 

When magnetic fields are close to equipartition, electrons should receive ap- 
proximately as much heat as ions, that is about O.lMc^. There are two different 
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mechanisms of electron heating. First, the viscous heat is evenly distributed be- 
tween electrons and ions (Gruzinov 1988, Quataert & Gruzinov 1999). Second, 
reconnection should be an efficient electron heater (Blandford 1998, Quataert &; 
Gruzinov 1999). 

If the electron heating is not negligible, as these papers suggest, the radiative 
efficiencies as small as 10~^ are impossible. Then the observed low luminosity 
of the non-AGN supermassive black holes should be explained by the small rate 
of accretion. Blandford & Begelman (1998) believe that the small accretion rate 
is explained by the wind loses. A different explanation is proposed in §3. 



3. The Accretion Rate 

The rate of spherical accretion onto an object of mass M is given by the Bondi 
formula 

^Bondi ^ PCsR\, (1) 

where p and Cg are the density and the speed of sound of the accreting gas at a 
large distance from the object, and Ra = 2GM/c^ is the accretion radius. For 
a supersonic wind, Cs is replaced by the wind velocity. The Bondi formula is 
widely used in astronomy. 

This formula was derived assuming that the flow is: (i) spherical, (ii) in- 
viscid, (iii) adiabatic. Although neither of these assumptions is realistic, the 
Bondi formula is a good order of magnitude estimate for the rate of accretion in 
many astronomically relevant cases (§3.1). We show, however, that in the case 
of low-density accretion onto a black hole, the Bondi formula can overestimate 
the rate of accretion by orders of magnitude (§3.2, §3.3) |^. 

For low-density quasi-spherical accretion onto a black hole, the accreting 
matter is a hot collisionless plasma. This plasma is an ideal heat conductor 
along the magnetic field lines. Heat conduction strongly violates the adiabaticity 
assumption used to derive the Bondi formula. 

In the absence of heat conduction, thermal energy of the plasma is lost 
through the event horizon. When heat conduction becomes important, some of 
this heat is transfered out, away from the hole. The temperature and pressure 
of the accretion flow increase over what have been the case for the adiabatic 
solution. The increased pressure opposes the gravitational attraction of the hole 
and impedes the accretion. 

In §3.1 we derive the Bondi formula and explain why it works for real (non- 
spherical, viscous, non-adiabatic, but relatively cool) flows. In §3.2 we show 
that if the effective adiabatic index of the flow 7eff = dlnp/dlnp is greater than 
5/3, the accretion slows down. We then show that heat conduction leads to 
7efT > 5/3, and derive a formula for the accretion rate in the presence of heat 
conduction. In §3.3 we present a more rigorous derivation of the same formula. 



^Most of the accretion disk papers treat M as an adjustable parameter. Here we actually 
calculate it, just like Bondi did, but under different assumptions. 
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When heat conduction is important, the rate of quasi-spherical accretion 
onto a black hole is 

M~MBondi (|^)°, (2) 

where Rs = 2GM/c^ is the Schwarzschild radius 0, and a ~ 1 is a model- 
dependent dimensionless number. 

To illustrate the astrophysical applications of the accretion rate formula (2), 
consider the case of Sgr A*, which is believed to be a 2.5 x 10® Mq black hole 
located at the center of our Galaxy (Genzel et al 1994, Ghez et al 1998). The 
accretion luminosity of Sgr A* is believed to be L ~ lO^^erg/s. When stellar 
winds in the vicinity of the galactic center collide and shock, they produce a 
gas of Cs ~ lOOOkm/s and p ~ IQ-^^g/cm^ (Coker & Meha 1997). The Bondi 
accretion rate of this gas is MRondi ~ 10^-^g/s. Thus MBondiC^ lO^^erg/s 
is five orders of magnitude higher than the actual luminosity. We (Gruzinov 
1998, Quataert & Gruzinov 1998) and others (Shvartsman 1971, Meszaros 1975, 
Blandford 1998) have argued that such low radiative efficiencies are unrealistic. 
It is then natural to assume that Sgr A* accretes at a much smaller rate (Bland- 
ford & Begelman 1998). If we use the estimate (2), and, for illustrative purposes, 
assume a = 0.4, the estimated accretion rate is reduced by a factor of 100. The 
radiative efficiency is then 0.1%. This radiative efficiency is still small, but it 
might be reasonable because electrons might be heated less than ions, and even 
hot electrons might radiate inefficiently when the plasma is rarefied (Quataert 
& Narayan 1998). 

3.1. The rate of spherical accretion according to Bondi 

The Bondi formula can be explained as follows. At a distance r ~ Ra from the 
massive object, the gravitational energy of the gas particles is comparable to the 
thermal energy. It is plausible that all these particles fall onto the black hole. 
Further assume that at r ~ Ra the radial velocity is ~ Cs and the density is 
~ p. We then get the mass accretion rate given by (1). 

Is the pressure build-up on the way into the massive object important? Can 
the pressure forces counter the gravitational attraction? We can answer this 
question by first neglecting the pressure forces at small radii, <C Ra- Then the 
radial velocity is equal to the free fall velocity, oc r~^/^. From mass continuity, 
the density is oc r~^/^. The pressure is oc r"^"^/^, where 7 is the adiabatic index 
of the gas. The acceleration from the pressure gradient is oc r~^^"'~^^^'^ . The 
gravitational acceleration is oc r~^. So long as 7 < 5/3, the pressure is not able 
to resist gravity. 

This qualitative derivation explains why the Bondi formula should work 
well if the accreting gas has an effective adiabatic index 7efr < 5/3. In this case, 
the pressure forces cannot impede accretion, and the gas from r ~ Ra must 
accrete at about the Bondi rate (1). Numerical simulations of a non-spherical 
accretion confirm the Bondi formula (Coker & Melia 1997). 



■^Relativistic effects are much smaller than the effects we are considering in this paper. We treat 
our black holes as Newtonian holes in space, and neglect the special-relativistic effects. 
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3.2. 7eff > 5/3 

Suppose that the accreting matter is a collisionless plasma. Suppose that ra- 
diation is neghgible. We will show that heat conduction increases the effective 
adiabatic index to 7eff > 5/3. Qualitatively, the heat conduction reduces the 
loss of thermal energy through the event horizon. The temperature of the flow 
increases faster than for an adiabatic plasma, which means that 7eff > 5/3. 
Quantitatively, TdS = dQ, or 

- pTvdrS = r-^dr{r'^KdrT). (3) 

Here p is the density, T is the temperature, v is the inflow (positive) velocity, s 

is the entropy per unit mass, n is the thermal conductivity. We assume that the 
thermal conductivity k is given by a "Shakura-Sunayev-like" formula 

K = m~^aprv, (4) 

where a ~ 1 is a dimcnsionless positive constant, m is the molecular mass. The 
Shakura-Sunaycv-likc prescription is plausible because the characteristic size and 
frequency of the turbulent structures should be ~ r and ^ v/r. The thermal 
conductivity might be suppressed for a supersonic inflow, but our accretion flow 
is subsonic. 

With s = m~^ln(T^/^/p), equation (3) gives 



pTv[~-^)+ar-\r'pvTy = 0, 



(5) 



where the prime denotes the r-derivative. For a stationary flow pvr^ = const, 
and 

'3T' p'' 



2T p 



+ a{rT')' = 0. (6) 



If a is small, we can use the adiabatic relationship in the last term, rT' = —T. 
Then equation (6) gives 

ri-" oc p, (7) 



that is 

5 -2a 
= 3^- 

Since 7eff > 5/3, pressure will impede the accretion. We now estimate the 
suppressed accretion rate. Assume that at the event horizon the temperature is 
close to m(P. The density at the hole ph is given by the adiabatic law 

1 ^ 



P Ws \Rs) 



)--\ (9) 



The rate of accretion is M ~ Rgcph, and from (8) and (9) we get (2). 
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3.3. A solvable model 

In this section we rigorously derive the accretion rate when heat conduction is 
important, equation (2). We assume a sphericahy symmetrical gas inflow. The 
radiative cooling is neglected. The flow is taken to be inviscid but non-adiabatic, 
i.e. we include thermal conduction and neglect viscosities. 

In terms of suppressing the rate of accretion, heat conduction is more im- 
portant than non- sphericity or viscosity. Non-spherical accretion is known to 
proceed at about the Bondi rate (Coker &; Mclia 1997). Viscosity suppresses the 
accretion rate, but only weakly, because when the inflow velocity decreases, the 
effects of viscosity decrease faster than the adiabatic effects (see also Turolla &: 
Nobih 1989). 

We treat the flow using a fluid approximation, but we should remember 
that in reality we are dealing with magnetized collisionless plasmas. However, 
the fluid approximation obeys the correct set of conservation laws, and should 
be close to reality. The qualitative conclusions of our analysis should be taken 
seriously. 

Consider a non-relativistic radial inflow of the 7 = 5/3 gas. The gas density 
p, the inflow (positive) velocity v, and temperature T depend on the radial 
coordinate r only, and satisfy the stationary equations of continuity, Euler's, 
and the thermal conduction. 

M = Airr^pv, (10) 
„„' = _W_^, (U) 

m \p 2 T J ^ ' ^ ' 

We will use dimensionless units with GM = 1, and the molecular mass m = 1. 

The boundary condition at the surface of the compact object is unknown, 
but it is irrelevant. For concreteness, assume T -|- /3rT' = at r = Rs, where (3 
is a dimensionless constant. At r = cxo, we assume T = 1, so that Ra ~ 1- The 
boundary condition for p is also irrelevant, because the system of equations is 
invariant under p \p, M AM. 

Our aim is to find the maximal possible value of the dimensionless mass 
accretion rate rh = M/p{oo), for which a smooth solution of the system exists. 
The ratio M/AI-Qondi is equal to this maximal value. Obviously m ~ 1 if Rs ~ 1. 
We need to find out how m depends on Rs when Rs <^ 1. 

Scale out M hy p ^ Mp/A-rr, insert (4) and (10) into (12): 

pvr"^ = 1, (13) 
vv' = -p-^ipT)' -r-'^, (14) 

T(^-l%)=a{rry. (15) 



P 2 T, 

Using (15), integrate (14) and obtain a system of two first-order differential 



equations 



6 



Andrei Gruzinov 



^(T-^^-^.^-r. (17) 



We need to find the dependence on Rs <C 1 of the maximal possible value of 
m = l/p{oo) for which a smooth solution of the system exists. For small r, 
we can neglect 5/2 in the right hand side of (16). Then, denoting r = rT, 
f = {r^p^)~^, X = logr, we obtain 

dr 

2a— = 2 - (5 - 2a)r - /, (18) 

dl ^ ^ 2 + 2dT/dx-5T-f 
dx T — f 

This system has a stable equilibrium point r = 2(5 — 2a)~-^, / = (it is assumed 
that a < 5/2). In the vicinity of this point, 

The system (18), (19) is an accurate approximation only for r ^ 1, i. e., for an 
x-duration 6x ~ log(l/i?5). The quantity / decreases according to (20) for an 
x-duration 6x = log{l/Rs) — const. Therefore the final density is proportional 
to (i?5)~°. Since rh = l/p{oo), the scaling law (2) is proven. To check the 
answer, we integrated equations (16), (17) numerically. 



4. The Density and Temperature of the CoUisionless Halo 

The standard Bondi solution and the new solution of §3 predict different densities 
near the black hole. In both cases the temperature is close to virial 

T{r) r^rupc'^^ (21) 
The density is given by the adiabatic. For the Bondi solution 

n{r)^n^(^^y\ (22) 

where noo is the density at a large distance from the hole (at accretion radius). 
When heat conduction is important, we use the effective value of the adiabatic 
index. This gives 

(T-> \ 3/2— a 
-f) ■ (23) 

For Sgr A*, under the assumption of §3, the new density is two orders of 
magnitude smaller. The density profile should, in principle, be observable. The 
deviations from the r^^/^ law will prove that collsionless halos around black 
holes are not ideal. 
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5. Discussion 



The low luminosity of the supermassive black holes existing in the centers of 
most galaxies (Richstone 1998) is traditionally explained by the extremely low 
radiative efficiency of the spherical accretion. We offer an alternative explana- 
tion: these black holes accrete at a rate much smaller than the Bondi rate. 

We have explained that electron heating and the radiative efficiency of the 
collisionless halos around black holes should be non- negligible. Then the ob- 
served low luminosity of many supermassive black holes should be due to the 
small accretion rate. 

We generalized the Bondi formula by including thermal conductivity. Ther- 
mal conductivity heats the flow, by transporting heat from the event horizon 
outward into the plasma. As the result of heat conduction, the pressure rises 
and impedes the accretion flow. 
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